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Abstract 

> 

Q\ • In this letter, the problem of optimal resource power allocation and relay selection for two way 

00 1 

T _j- ■ relaying cognitive radio networks using half duplex Decode and Forward (DF) and Amplify and Forward 

(AF) systems are investigated. The primary and secondary networks are assumed to access the spectrum 

m ! 

at the same time, so that the interference introduced to the primary network caused by the secondary 
network should be below a certain interference threshold. In addition, a selection strategy between the 
AF and DF schemes is applied depending on the achieved secondary sum rate without affecting the 
quality of service of the primary network. A suboptimal approach based on a genetic algorithm is also 
presented to solve our problem. Selected simulation results show that the proposed suboptimal algorithm 
offers a performance close to the performance of the optimal solution with a considerable complexity 
saving. 

Index Terms 

Cognitive radio network, two way relaying, relay selection, genetic algorithm. 



I. Introduction 

Cognitive Radio (CR) has recently attracted enormous attention in wireless communication 
networks [1]. It is considered as a promising solution towards a more efficient usage of the 



radio spectrum. The idea of CR spectrum sharing is to allow unlicensed users known also as 
Secondary Users (SUs) to utilize the spectrum band allocated by licensed users known also as 
Primary Users (PUs) at the same time. In order to protect the PUs, the interference due to the 
SUs should be kept under a certain interference temperature limit. 

On another front, there has been recently a great deal of interest in two way relaying networks 
[2], [3]. The transmission process in this relaying technique takes place in two time slots. In 
the first slot, the source and the destination transmit their signals simultaneously to the relay. 
Subsequentaly, in the second slot, the relay broadcasts its signal to the terminals. Two widely 
relay protocols are used in practice, the namely Amplify and Forward (AF) protocol, which 
amplifies the received signal first, then broadcast it to the destination, and the Decode and 
Forward (DF) protocol, which decodes the received signal to remove the noise before transmitting 
a clean copy of the original signal to the destination. For instance, the work presented by Chen 
et. al in [2] deals with multi access two way relaying network case, while in [3] the authors 
show analytically and via simulation that two way relaying outperforms one way relaying in 
terms of energy efficiency. Furthermore, the relay selection and power allocation problems for 
AF protocol in cooperative one way and two way relaying CR have been investigated in [4] and 
[5], respectively. The best relay selection in two way relaying depends on two factors, end to end 
channel conditions, and the presence of the Primary Network (PN) according to the interference 
constraints . On the other hand, prior work in the literature has studied adaptive relaying which 
allows the switching between AF and DF protocols depending on the Signal-to-Noise-Ratio 
(SNR) [6]. However, to the best knowledge of the authors, the relay selection problem in two 
way relaying CR networks using DF protocol has not been discussed so far as it is the case for 
the AF protocol. 

In this letter, a best relay selection scheme for two way relaying CR with half duplex case 
and channel reciprocity is considered. In the AF protocol, the relay broadcasts the amplified 
copy of the received signal to the terminals, i.e., the noise gets amplified too. On the other hand, 
in the DF protocol, the relay regenerates clean signals from the received signal and transmits 
the re-encoded message to the terminals. More specifically, the main contributions for our new 
proposed scheme can be summarized as follows: 

• Formulate a new relay selection scheme in two way relaying CR system which selects 
between the DF and AF protocols depending on the higher Sum Rate (SR) achieved by the 



Secondary Network (SN) without affecting the Quality of Service (QoS) of the PN. For 
that reason, additional interference constraints are considered in the optimization problem 
for both time slots (it is assumed that I th is the same in each time slot). 

• Derivation of the optimal transmits power and relay power that maximize the cognitive SR 
of the system. 

• Using dual decomposition and subgradient methods for both AF and DF techniques in order 
to solve the SR maximization problem and select the best relay with the best technique. 

• Design a practical low complexity suboptimal approach based on Genetic Algorithm (GA) 
to solve the formulated optimization problem [7], and compare it with the optimal and 
Exhaustive Search (ES) solutions. 

Generally, in one way relaying, it is assumed that at high SNR the relay can decode perfectly, 
so it achieves higher SR using the DF protocol. On the contrary, for low SNR the higher SR 
can be achieved using the AF protocol. However, the results provided in Section V show that 
in two way relaying at high SNR the DF protocol becomes as a bottleneck in the first phase, so 
higher SR can be achieved using the AF protocol. On the other hand, for low SNR, the relay 
with the DF protocol achieves higher SR. 

The rest of this letter is organized as follows. Section II gives the system model. The problem 
formulation and the optimal algorithm are described in Section III. The suboptimal scheme 
is presented in Section IV. Simulations and numerical results are demonstrated in Section V. 
Finally, the letter is concluded in Section VI. 

II. System Model 

In this section, the best relay selection problem for CR two way relaying is investigated. The 
SN is constituted of a cognitive Mobile User (MU), a Cognitive Base station (CB), and M 
Relay Stations (RSs) as illustrated in Fig.l. It is assumed that there is no direct link between 
the cognitive terminals and the single relay principle is applied to select the best relay. During 
the first time slot known also as the Multiple Access Channel (MAC) phase, the CB transmits 
its signal to RSs with power denoted Pqb- Concurrently the secondary MU transmits its signal 
to RSs with power denoted P s . This causes two interferences to the PU. In the second time 
slot known also as the Broadcast Channel phase (BC), the selected RS broadcast its signal. This 
phase also causes interference to the PU from the active RS. 




Fig. 1. System model of the cooperative two way relaying cognitive radio system. 



We assume that all the channel gains are perfectly known at the communication nodes. All 
channel gains for the network can be adopted by assuming channel reciprocity and classical 
channel estimation approaches [8]. The interference between the PN and SN is studied in 
Section III. Also, we assume that the PN and SN access the spectrum at the same time. Without 
loss of generality, all the noise variances are assumed to be equal to <x^. Finally, selection strategy 
between the AF and the DF protocols is applied in order to achieve the maximum SR of the 
SN without affecting the QoS of the PU measured by I th . 

III. Problem Formulation 

Let us define Rdf and Raf as the achievable secondary SR for the DF protocol, and the 
achievable secondary SR for the AF protocol, respectively. 

The Optimization Probleml (OP1) for a single relay selection can be formulated as 1 

(OP1): m* = argmax max pRdf + (1 — p)Raf, (1) 

me{l:M} R 

s.t 0<P s <Ps, (2) 
< Pcb < Pcb, (3) 



'For simplicity and uniformity we use the mathematical notations depicted in Table I. 



0<P R ,m<P R , Vm=l,...,M, 



(4) 



interference constraint in first time slot 



hP S + f±PcB 



<Ith, 



(5) 



interference constraint in second time slot 



fcP R , m <Ith, Vm = l,..,M, (6) 

where P$, Pcb, and Pr, are the peak transmit power of the secondary MU, CB, and m-th RS, 
respectively. In (1), p is a constant equal to either zero for the AF protocol or one for the DF 
protocol, and the channels coefficients are given in Table I. Let x\ and x 2 are the symbols 
transmitted by the MU and CB respectively. It is assumed that E(|xi| 2 ) = E(|x 2 | 2 ) = l 2 . In the 
first time slot, the received signal at the m-th relay is given by 

r m = ^fP S h i i- R) x 1 + v^W4 CB ~ R) ^ + z m , (7) 
where z rn is the additive Gaussian noise at the m-th relay. 



TABLE I: Symbol Notations 



Symbol 


Notation 


Complex Channel Gain between 


h 


\hlS B - R) \ 2 


CB and RS m 


h 


\hiS- R) \ 2 


MU and RS m 


h 


| fe (S-P)|2 


MU and PU 


h 




CB and PU 


h 


|^- P) | 2 


RS m and PU 



In order to simplify the formulated OP1, we solve it time slot per time slot. During the BC 
phase, the power allocation at the m-th relay depends essentially on two constraints: the peak 
power constraint (4) and the interference constraint (6). For this reason, the optimal relay power 
can be expressed as 

Pl m = mm(pJ-fY Vm = l,..,M. (8) 



2 E(-) denotes the expectation operator. 



The optimization problem during the MAC phase is therefore given by 

(OP2): m* = argmax max pR DF + (1 — p)R AF , (9) 

me{l:M} R 

s.t (2), (3), (5) (10) 

We can decompose the OP2 outlined above into parallel subproblems using single relay principle, 
i.e., each independently solvable for a different relay, then we select the relay that achieves 
maximum SR. 

The dual subproblem associated with OP2 can be written as [9] 

mm g(X), (11) 

where A is a lagrangian vector contains the Lagrangian multipliers in the system. The dual 
function g(X) is defined as follows 

g(\)= max £{\, P s , P CB ). (12) 

P,S>0,PCB>0 

A. Amplify and Forward Protocol 

In this protocol, The relay amplifies the received signal by w m , then the received signal at 
the terminals can be expressed as 

r m ,s = w m h^~ R) r m + z CB , 

(13) 

r m ,CB = W m ll { m B R)r m + Z S , 

where z C b and zs are the additive Gaussian noise at the terminals. By using the perfect 
knowledge of the channel gains and channel reciprocity, the terminals can remove the self 
interference by eliminating their own signals. Thus, the SNR at MU and CB are given by 

PcB\ \ 2 f2h Ps\w m \ 2 f 2 fi .... 

<(Fm| 2 /2 + 1) 0*{\w m \ 2 h + 1) 

respectively. The relay power of the m-th relay node can be expressed as 

P R)m = E(\w m r m \ 2 ) = (P5/2 + PcbIi + °l)\wm\ 2 . (15) 

By substituting the value of \w m \ 2 from (15) into (14), the SNRs become 

PcBP% m hh 



Jm,S 



lm,CB 



* 2 n(P* R , m fl + Psf2 + PcBfl+<T 2 n y 
_ PsPR^mfzfl 

= <(PR, m f2 + Ps.f2 + PcBh + a 2 n ) 



(16) 



The achieved SR for AF protocol of two way relaying can be written as 

RAF = ^ l0g 2 (l + 7m,s) + \ l0g 2 (l + 7m,CB)- (17) 

Due to the non-convexity of the formula in AF protocol, a convex approximation when the 
system operates at high SNR region is presented [5]: 

R AF ~ i log 2 (7 m , 5 ) + ^ log 2 (7 m , C B). (18) 

When p = and due the fact that the logarithmic function is a monotonically increasing function 
of its arguments, the OP2 is equivalent to the following 

(OP2): m* = argmax min , (19) 

me{l:M} lm,CB-lm,CB 

s.t (2), (3), (5) (20) 
Thus, the Lagrangian of AF protocol can be written as 

£af = As(Ps — P s ) — Xcb{Pcb — Pcb) — ^i{hPs + IaPcb — hh), (21) 

lm,CB-lm,CB 

where As, \cb, and Ai represent the Lagrangian multipliers related to the peak power at the 
source, peak power at the destination, and interference constraint in the first time slot, respec- 
tively. By applying the Karush-Kuhn-Tucker (KKT) optimality conditions [9], we obtain 

^ AF n A ^AF n 

—— = and — — = 0. (22) 

rs rCB 

Direct calculation yields 



P * CB kI Ua 2 + (acb +*u)p s pg, m ftn) (24) 

where A = p£j 2 f\ + P^ m (f 2 a^ + j\a 2 n ) + I% B fi + P CB {V^l) + ^U^WA/i + f!) + o£, 
^ = ^ m / 2 /i + I% m {M + Vn) + ^l/l + ^5(2/2^) + P* R , m Ps(f2fi + /I) + < and (x)+ 
denotes a maximum between x and zero. 



B. Decode and Forward Protocol 

Prior works in the literature have studied the sum rate for two way relaying with DF protocol 
[10]— [12]. The max SR of the DF protocol can be expressed as 

1 



Rnp = - min 
2 



min{ J R 1 , R 3 } + mm{R 2 , R 4 }, R 5 



(25) 



where R t = log 2 ^l + ^ffi^j, R2 = log 2 (l + ^js^ j, denote the rate from the source and 
the destination to the relay in the first time slot, respectively, R 3 = log 2 ( 1 + ^p- ) , Ra = 
log 2 ( 1 + J , denote the rate from the relay to the source and to the destination in the second 
time slot, respectively, and R 5 = log 2 ^1 + p cBh+Psh ^ d eno t e s the max SR can be achieved in 
both time slots. 

It is assumed that the relay node decodes the high SNR signal (Down-Link (DL) signal) first, 
then decodes the other signal (Up-Link (UL) signal) after subtracting the decoded signal. For 
this reason additional Lagrangian multipliers are considered for UL and DL. When p = 1, the 
Lagrangian of OP2 can be written as 

£df = (1 - K - 1 + A d )|log 2 (l + - A d )±log 2 (l + p cBfi+Psf2 ) _ Xs{ p s _ p s) _ 

^cb(Pcb — Pcb) — ^i{hPs + UPcb — hh)- 

(26) 

where \ u and \ d are the dual variables associated with the UL and DL rate constraints, respec- 
tively. Letting a = 2. In 2 and applying the KKT optimality conditions, we obtain 

= and = 0. (27) 

Ps Pcb 

Direct calculation yields 

P ° B ~ W^U + Xcb) ~ (28) 

Kl Pf + k 2 P* + k 3 = 0, (29) 

where Kl = (X 1 f 3 + \ s )fln 2 = f 2 (2a 2 n + P CB /i)(Ai/ 3 + \ S ) ~ and k 3 = (1 - 

\ d )PcBh% + K(Ai/ 3 + A s ) - (1 - \u)^)(a 2 n + Pcsfi). By substituting (28) into (29) and 
after simplification, we obtain the optimal source power as the following 

( \ 1 \ \( hflO-- x d) \ „2/i x ,\( ( A l/3 + A s)/l /2-v 
P* = I [ d ~ uK ° 2 ^lf4 + *CB)> n( dK ^lU + X CB)--> 1 HQ) 
S I (l-A d )/2/l(Ai/ 3 +A s ) | (A d -A u )/f-(l-A„)/ 2 I 



a(A 1 / 4 +A cs ) 



C. Dual Problem Solution 



The dual problem of OP2 can be solved by using the subgradient method [13]. Therefore, to 
obtain the solution, we can start with any initial values for the different Lagrangian multipliers 
and evaluate the optimal powers, then update the Lagrangian multipliers at the next iteration as 



\t+i _ \t 



5(t) 



p* 



(31) 



\t+i 

A CB 



A CB 



-5(t) 



PCB — P> 



CB ) 



(32) 



AJ+ 1 = A* - 5(t) [l th - [hP* s + UP* CB ) 



K +l = a* - s(t) [I i og2 (i + %A) _ i log2 (i + 3£) 



(33) 



(34) 



P*sf2 + P* CB h' 



I log 2 ( 1 - 



(35) 



where 5(t) is the step size updated according to the nonsummable diminishing step lengths 
policy [13]. Using the subgradient method, the updated values of the optimal powers and the 
Lagrangian multipliers are repeated until convergence. The implementation procedures to solve 
the OP2 is described in Algorithm 1. 

IV. Suboptimal Algorithm 

The optimal solution for our non linear OP2 sometimes is difficult to solve due to its high 
computational complexity. Therefore, in order to solve the problem efficiently, we propose a low 
complexity suboptimal approach in discrete domain to find suboptimal solution. In the MAC 
phase, we need to find the optimal power allocation over the terminals (i.e., P s and Pqb ) in 
order to maximize the SR of SN without interfering with the PUs. 

In this section, we propose a heuristic GA with discrete number of power levels from zero 
to the peak power budget. In fact, each terminal can transmit its signals using one of the 
power levels between and peak power budget, i.e., [p s e jo, j^L, ^ N ^[ s , ^s})' 



Algorithm 1 Optimal Power Allocation and Relay Selection 



- Input: I th , P s ,Pcb,Pr, M, f lt f 2 , h, h, h- 

~ Rmax = 0' 

for m — 1 : M do 

- Initialize the Lagrangian multipliers A, Pcb, and p = 0. 
while p = do 

- Solve problem (23) to obtain P§, P s = P|. 

- Solve problem (24) to obtain P£ s , Pcb = P^b- 

- Update A using subgradient method based on (31) - (33). 

- Until Required precision is satisfied or reach maximum iteration, 
end while 

- Find Raf using (18) 

- Initialize the Lagrangian multipliers A, and p=l. 
while p = 1 do 

- Solve problem (30) to obtain P*, P s = P*. 

- Solve problem (28) to obtain P£ B . 

- Update A using subgradient method based on (31) - (35). 

- Until Required precision is satisfied or reach maximum iteration, 
end while 

- Find Rdf using (25) 

- Rmax — max(P J 4j?, Rdf) 
end for 

- Find m* s.t R opt = max R max 



and [Pcb e |0, jf^, , (JV ~y B , p CBjj where N is the number of quantization levels. 
In this way, the transmitters have more flexibility to allocate their powers in the case where 
continuous power distribution is not available. The GA tries to find the optimal binary string 
that maximizes the SR expressed in (9). At the beginning, we represent the discrete quantization 
values of P s and P C b as iV binary strings each of length K, where 3 K = |~log 2 (iV)]. The 
binary representation set of P s and Pcb are denoted as Ss and Scb, respectively. The algorithm 



3 [a;] denotes the smallest integer not less than x. 



concatenates S s with S C b to produce an initial population set S of N elements and each with 
2K bits, where the first K bits represent the equivalent binary string for P s and the last K 
bits represent the equivalent binary string for Pcb- Initially, the GA computes the SR of all 
elements in .So using (9), then maintain the best y strings e S to the next population and from 
them, generate y new strings by applying crossovers technique to form a new population S. This 
procedure is repeated until reaching convergence (i.e., SR remains constant for several successive 
iterations) or until reaching the maximum generation number /. Details of the proposed GA are 
given in Algorithm 2. 

Algorithm 2 Proposed Genetic Algorithm 

- Input: I th , P s ,Pcb,Pr, M, f u f 2 , h,h, hJ- 

Rrnax 0- 

for m = 1 : M do 

-P^ m = mm (Pr,%). 

- i = 1, Ri = 0, and generate an initial population set S. 
while (i < I or not converge) do 

for n = 1 : N do 

if (interference constraint is satisfied) then 

- = Compute the sum rate using (9). 
else 

- i?W = 0. 
end if 

end for 

- Rf = max( J R). 

- Maintain the best ^ strings G S to the next population and from them, generate y new strings by 
applying crossovers to form a new population S. 

- i = i + 1; 
end while 

- Rimax = max( J R 7 ). 
end for 

- Find to* s.t R opt = max R ma x- 



The formulated OP2 can be, of course, solved via an ES algorithm by investigating all 
possible combinations of the transmitters power and select the best combinations that satisfied 



2 

the interference constraint This algorithm requires M { 2 )(N ~~ 1)* — 0(MN 2 ) operations 

i=o 1 

[14]. However, our proposed GA requires MNI operations to reach a suboptimal solution. 

In the proposed algorithm, our goal is to maximize the SR of the SN without interfering with 
the PU. The last step in our proposed algorithm is selecting between the AF and DF protocols 
depends on the higher achieved SR. Hence, our proposed algorithm is able to reach a suboptimal 
solution with a considerable complexity saving. In addition to that, simulation results in Section 
V show that by increasing N, our proposed GA achieves almost the same performance as the 
optimal solution. 

V. Selected Simulation Results 

In this section, some selected simulation results are performed to show the benefits of our 
system. We assume a single cell subject to a small scale Rayleigh fading, consisting of one PU 
and a SN constituted by one CB, one secondary MU, and M = 4 relays. The variance o 2 n is 
assumed to be equal to 1CT 4 . We also assume that the transmit power constraint of MU, CB, 
and each RS are equal to P bar . The proposed GA is applied under the following settings: the 
crossover point is chosen randomly between 1 and 2K for each binary string, and we run the 
GA at most 10 times. 

The advantage of adaptive relaying strategy is depicted in Fig. 2. The adaptive strategy can 
switch between the DF and AF protocols according to the best performance. It is worth men- 
tioning that, in the high SNR regime, adaptive relaying uses the AF protocol. On the other hand, 
for the low SNR region, adaptive relaying uses the DF protocol. Fig. 2(a) plots the SR versus 
peak power P bar , while Fig. 2(b) plots the SR versus interference threshold, for different values 
of I t h and MU peak power, respectively. In general, the results suggest the usage of the AF 
protocol when both P bar and I t h are large. This can be justified by noticing that the SR value 
of the DF protocol becomes as a bottleneck for the first phase in the high SNR regime. 

Fig. 3 shows a comparison between the performance of the proposed GA with the optimal and 
ES solutions. We plot the achieved secondary SR versus P bar for different values of I th = {20, 5} 
dBm and different relaying protocols. We can notice that, in the low P bar region, the proposed 
GA, the optimal solution, and the ES have almost the same sum rate, while in the high P bar 
region, a gap between these methods is observed. This gap is increasing with higher P bar . values. 
This is justified by the fact that starting from a certain value of P bar the GA can not supply the 



selected relay with the whole power budget. Hence, the selected relay transmits its signal with 
one of the quantized power levels. In fact, with high values of P bar , the constraint (5) can be 
affected. For this reason, we introduce the discretization set to get more degrees of freedom by 
increasing TV and as such enhance the SR. For instance, Fig. 3(a) and Fig. 3(b) plot the secondary 
SR for I th = 20 dBm for DF protocol and AF protocol, respectively. It is shown that the GA 
achieves almost the same SR reached by the optimal solution. However, when I th is reduced, 
we notice a degradation of the GA performance at large values of P ba r as shown in Fig. 3(c) and 
Fig.3(d). 

The same interpretation is applied on Fig.4 in which the achieved secondary SR is plotted 
versus the interference threshold for both relaying protocols. In this figure, for fixed Pb ar the 
performance of the GA is close to the optimal solution for large I th . One can see that, a gap 
between the methods is noticed in the low I t h region. This can be justified by the fact that in 
this region the GA can not reach the maximum power budget due the small value of I t h- Hence, 
the GA tries to transmit with one of the quantized power levels. However, It can be shown that 
when N — > oo, the proposed GA achieves the performance of the optimal solution. 

VI. Conclusion 

This letter presented an optimal power allocation and a relay selection scheme for two way 
relaying cognitive radio networks using the dual decomposition method. The idea of this scheme 
is to maximize the SR of the cognitive network taking into consideration protecting the PUs 
from the interference caused by the SN. Due to the high computational complexity of the optimal 
solution, a suboptimal heuristic algorithm is presented. The suboptimal solution based on the 
GA is able to achieve the same performance of both the ES and optimal solutions with a much 
less complexity. Furthermore, the performance of the DF and AF schemes, and the impact of the 
power and interference constraints are illustrated for different interference thresholds and peak 
power constraints. Finally, the advantage of the adaptive relaying protocol is shown by switching 
between the DF and AF protocols. 
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Fig. 3. The achieved SR of the proposed GA, the ES algorithm, and the optimal solution with different values of I t h, and N 
versus Ptar, for (a,c) DF protocol, (b,d) AF protocol. 
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Fig. 4. The achieved SR of the proposed GA, the ES algorithm, and the optimal solution with different values of Pt, ar , and iV 
versus I t h, for (a,c) DF protocol, (b,d) AF protocol. 



